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Abstract. We study numerically a two-dimensional version of the Burridge-Knopoff 
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1 Introduction 



A mathematical model of an earthquake fault, originally introduced by Burridge 
and Knopoff [|T|, has been studied extensively in recent years, see e.g. [| U [5], as 
an example of a self-organized critical system |], [?|. The motivation has also been 
to understand better the mechanism of earthquake generation as well as the search 
for methods to predict earthquakes |8], ||. The quantities of main interest have been 
the scaling law (Gutenberg-Richter law |H|) relating the frequency of earthquakes 
to their size, other scaling laws as well as spatial and temporal correlations. 
In 
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the one-dimensional model is studied in detail by integrating numeri- 
cally the equations of motion. In [12j a simplified model with qualitatively similar 
behaviour but simpler time development, corresponding to coupled maps, was intro- 
duced. This model, which is analogous to sandpile models J0|, was studied in more 
detail in 0, [14] . Some special two-dimensional models are studied in ||, [15], [16| , see 
|T7| , [T8| | for older work in the same vein. The purpose of this article is to generalize 
the study of |14[ to two dimensions and compare the one and the two dimensional 
cases. 

We begin by defining the model and the observables of interest. Then we describe 
the results of extensive simulations of the size distribution, spatial and temporal cor- 
relations as well as the structure of the slip zone. We briefly describe the algorithm, 
which is essentially the one of |19|, and finally we compare our results with previous 
work and discuss some outstanding questions. 



2 The model 

We describe here the two-dimensional version of the Burridge-Knopoff model and 
it will be clear from the discussion how to generalize it to higher dimensions. The 
higher dimensional analogue does of course not have any straightforward physical 
interpretation. We restrict ourselves to the homogeneous model. The model can be 
described as a two-dimensional regular array of blocks of equal mass, which slide, 
subject to a velocity dependent friction, on a plane which can be thought of as one 
side of an earthquake fault. The blocks are labelled by pairs of integers (i, j), where 
i = 1, . . . Ni and j = 1, . . . , N 2 , such that x it j is the coordinate of the block labelled 
by {hj)- Each block is connected to its nearest neighbours in the first coordinate 
direction by a spring with spring constant k and to its nearest neighbours in the 
second coordinate direction by leaf springs with a spring constant k! . We only 
consider the isotropic case where k' = k. In addition, each block is connected by 
a leaf spring, with spring constant k", to a sliding plane which lies parallel to the 
plane of the blocks and moves with velocity v p . This plane can be thought of as the 



other side of the fault. For a discussion of boundary conditions we refer to |16|. For 



convenience we choose units such that k" = 1 and v p = 1. 

If we start the system at rest in a configuration where the elastic force on each 
block is smaller than the static friction force, then the blocks remain at rest until 
the elastic force on one of the blocks exceeds a threshold value f t h which equals the 
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static friction. Then this block starts to slide, with increasing velocity to begin with, 
since the friction is a decreasing function of velocity for small velocities in physically 
realistic models, see [[HJ. The motion of the first block may cause the force on some 



of the neighbouring blocks to exceed threshold so they begin sliding and so on. In 
some cases the motion of the first block unleashes an avalanche of sliding blocks and 
this is regarded as a large earthquake. Eventually all the blocks come to rest again, 
at least if v p is small enough compared to the natural velocities in the spring-block 
system. 

We normalize the coordinates such that Xij is the displacement of the block 
from its equilibrium position at time t — 0. The total nonfrictional force fij on this 
block at time t due to the springs is then given by 

fi,j ip^i>j t) ~\- k ^ ] (^n,m ("0 

(n,m) 

where the sum runs over the four nearest neighbours of the site We do not 

study the Newtonian equations of motion for this system (see 0) but rather use 
the simpler discrete dynamics ([12, ||) which we now describe. The elastic force 



increases linearly with time on each block, so there is generically a unique first block 
which reaches the threshold for moving. We let this block move along the first 
coordinate axis (the direction of the driving force) until the elastic force on it has 
reached the value 

tij = oil,,). (2) 

where is a decreasing function which describes the force relaxation in the system. 
This motion is assumed to take place in zero time and meanwhile all the other 
blocks remain at rest. For earthquake dynamics this is an excellent approximation. 
The slip of the block gives rise to an increase in the force on its four nearest 
neighbours by an amount ~A(/jj — 4>(fi t j)), where 

' ^ (3) 



l + 4fc 



For a (i-dimensional system 4 is replaced by 2d. In general this force relaxation is 
modified at the boundary. The parameter A is a measure of the stiffness of the 
system or, if one prefers, the amount of conservation of elastic stress in the system, 
A = 1 corresponding to a conservative system. 
In practice we take so that 

<f>(fth) = fth - Sf, (4) 

where Sf > is the smallest possible force relaxation. By making one more choice 
of units we may assume that f t h = 1- Using the linear relationship between the 
coordinates of the blocks and the elastic forces we can easily compute the slips of 
blocks during earthquakes. 

If the force relaxation function 4>(f) decreases fast for / close to fth, then a small 
slip is easily amplified as we discuss in more detail below. In our simulations we 
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follow & [12] and use = ip ai where 



ip a (x) 



a(x-f th ) + (2f th -5f) 



-fi 



tii 



(5) 



and a > 0. Note that ip' a (fth) = —ol so the size of a is a measure of how easily 
small slips can trigger large ones. We also discuss briefly the case = considered 
in B Il6| and in more detail in ]19 . 



3 Observables 

We now describe the quantities that have been of main interest in the study of the 
model. In addition we introduce a natural notion of a Hausdorff dimension for the 
slip zone. 

If Xij and x\ 3 - are the coordinates of the blocks before and after an earthquake, 
its moment is defined as 

Ni N 2 

^ = EEKr^ ( 6 ) 

t=l j=l 

There are two other quantities of interest in describing the size of an earthquake, 
the number of participating blocks which we shall denote by N and the energy E 
realeased which is obtained by comparing the potential energy in the springs before 
and after an event. 

The magnitude of an earthquake is defined as the logarithm of its "size" so 
different notions of "size" give rise to different notions of magnitude. In particular, 
we shall see that the moment magnitude \i = logM, does not obey a Gutenberg- 
Richter law in two dimensions except for very special choices of the parameters. 

Before proceeding let us discuss qualitatively the main features of the events 
observed in the case <f> = ip a . For a soft system the blocks move independently 
and only make the minimal jumps corresponding to the force relaxation Sf. In 
fact, the system is periodic when it is sufficiently soft; if a given block moves at a 
certain time, then all the other blocks move before the given block moves again. 
At a sharply defined value of A = A events with more than one block begin to 
occur in the system and there is no upper bound on the size of events that can 
occur, see Fig. 1, where this scenario is illustrated for the one-dimensional system. 
The precise value of Ao depends on a and the transition is sharper the larger the 
value of a. In fact our data are consistent with a discontinous change in the size 
of the largest events for a > 3. We interpret this transition as being due to the 
onset of instability in the periodic state. Once this state becomes unstable and the 
motion of a block can trigger a motion of its neighbours, there is little that prevents 
earthquakes from spreading over the entire system. As we increase A further the 
dynamics becomes more chaotic f2~lf , the state of the blocks more irregular and 
earthquake fronts are more likely to get trapped in "valleys" in the spring-block 
system. As we increase A still further the size of the maximal events passes through 
a minimum and eventually the system becomes so stiff that almost any motion of a 
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block can trigger a systemwide slip. In Fig. 2 we illustrate the same phenomenon 
for the two-dimensional system. 

Let us discuss this in more detail. The displacement of the first block to move 
in an earthquake is always 

0i = Sf(l - A). (7) 

Assuming that all blocks in an infinite system are at threshold the displacement of 
the i-ih block in an earthquake is inductively given by 

«, = l-A + ^« M -(l-A)*(l + ^ij). (8) 



This result was first derived for the one-dimensional system in |L4[] but is in fact 
independent of the dimensionality of the system. If is a decreasing function, then 
the sequence 6i is easily seen to be increasing and bounded so it converges to a value 
8* which is the maximal slip of a single block in an earthquake. 

Let us study the stability of the periodic state in a soft system more closely. For 
simplicity we discuss the one- dimensional case but the mechanism is the same in 
two and higher dimensions. We assume that the system is in a state where all events 
are minimal. Then we have for all i, 

1 - Sf < fi < 1, (9) 

where fi is the total force from the springs on block %. Let us assume that block 1 
is going to move next and fx = 1. We now also put f^ to threshold and we want to 
derive a sufficient condition on a in ip a and A so that block 3 will also move. We 
have 

fi^l-Sf (10) 

and 

h » K = i + \*Sf. (ii) 

Since 

1 - Sf < f 3 < 1 (12) 
a sufficient condition for the movement of the 3rd block is 

^A(l-0(l + iA5/))>5/. (13) 



With <f) = ip a this condition is equivalent to 

(2-5ff 4 



+ 2<aA. (14) 



Sf(l - f ) Sf 

Since 8f <C 1 it follows that the above inequality is roughly equivalent to 

4 f^- - l) +2 < aA (15) 
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which agrees quite well with the location of the onset of the instability in the qA 
plane, cf. Fig. I. 

We now introduce the correlation functions. Let F be a function of time and 
the coordinates Xij of the blocks. We refer to such functions as observables. The 
average < F > of an observable F is defined as 

<F>= lim I [ T Fdt, (16) 

T^oo L JO 

provided the limit exists. In the simulations we choose of course a finite value of T 
which must be taken sufficiently large, depending on the observable under study as 
well as the system size. 
We define 

€i,j(t) = Xij(t) - v p t. (17) 

The average < > is well defined since the long time average of x it j is proportional 
to t with coefficient v p . 

The spatial correlation function G(i) is defined as 

G(i) =< £i+jVi/2,jV 2 /2£iVi/2,Ar2/2 > — < £i+JVi/2,Ar 2 /2 >< (,N 1 /2,N 2 /2 > ■ (18) 

Here we have assumed that Ni and N 2 are even numbers and chosen the center of 
the system as a reference point in order to minimize finite size effects in the case 
of free or open boundary conditions. In the case of periodic boundary conditions it 
does of course not matter which subtraction point we choose. A more general spatial 
correlation function is obtained by separating the two £ variables by an arbitrary 
lattice vector. This more general correlation respects the symmetry of the lattice 
but does not reveal any features different from those of G(i). 

We define the time correlation function (or autocorrelation function) C(t) as 

C(t) = hm I [ T Z(t> + t)C(t') dt>- < i > 2 , (19) 

where £ = t, Nl /2,N 2 /2- Clearly C(0) = G(0). 

It is of interest to consider the shape of the slip zone in earthquakes. For this 
purpose we introduce the notion of the Hausdorff dimension S of the slip zone. It is 
defined as 

, IniV 

S = ^D < 20) 
where iV is the number of blocks that participate in the event and D is the length of 
the diagonal of the smallest rectangular box, with sides along the coordinate axes, 
which contains all the blocks that participated in the earthquake. This definition 
is by no means unique but we find it convenient to work with. All reasonable 
definitions should be equivalent in the limit of large iV and D. 



4 Numerical results 

In this section we describe the results of simulations for the relaxation function 
■03- Most of the results are for a 200 x 200 block system with periodic boundary 
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conditions. We begin by considering the size distribution of earthquakes. We let 
R(fi) denote the density of events per unit time with moment magnitude fi. There 
are three regimes in this distribution for all values of A and a. We discuss only 
the case of a = 3. Other values of a in the range from 2 to 5 do not lead to a 
qualitatively diferent picture. 

For small events the discretization shows clearly up. In the intermediate range 
from fi pa —2 to /i ~ the distribution falls off but not linearly except for A pa 0.82. 
For smaller A the distribution is concave in this region while it is convex for larger 
values of A, see Figs. 3-5. This is quite different from the one-dimensional model 
where the size distribution is linear to a good approximation in the intermediate 



range ||14|| . The distribution of earthquakes as a function of the energy is qualitatively 



similar. However, the distribution as a function of the number of participating 
blocks is linear to a very good approximation for A > 0.8 with a slope which varies 
continuously with A, see Figs. 6-7. For smaller A the distribution becomes concave. 

For A < 0.65 almost all events are minimal but around A = 0.65 the instability 
discussed in the previous section sets in, cf. Fig. 2, and for A = 0.67 some events 
extend over the entire system with almost no events of an intermediate size. As A 
increases further events of intermediate size begin to appear and finally for the very 
stiff system the peak at large events begins to rise again. 

The behaviour of the size distribution is also reflected in the correlation func- 
tions. Let us first discuss the spatial correlation G(i), see Fig. 8. For small A we 
find a random function of small amplitude since the system is periodic. As events 
involving more than one block begin to occur we find a sharply peaked function 
at which explodes at the "phase transition" point A = 0.67 to a broad function 
whose decay is not very well determined due to finite size effects and the very large 
fluctuations in the simulation at this value of A. For larger values of A the functions 
narrows until the prevalence of large earthquakes reaches a minimum at A pa 0.82. 
It broadens again as A is increased further. In an infinite system the correlation 
should decay monotonically to zero. We see marked finite size effects for A = 0.67 
as well as for A = 0.95. This is in a qualitative agreement with the behaviour of the 
spatial correlation in the one-dimensional model and the value of G(0) is roughly 
proportional to the maximum displacement of blocks in a single earthquake, cf. Jl4| . 

The autocorrelation function C(t) is considerably more difficult to calculate in 
two dimensions than in the one-dimensional case. As far as we can see it has 
qualitatively the same behaviour in one and two dimensions as long as the system 
is not too stiff. The correlation decays rapidly with well defined oscillations whose 
frequency is well explained by the maximum slip of a single block and the driving 
velocity, see Fig. 9. The data does not allow us to say anything about the decay of 
the envelope. For a stiff system (A pa 0.95) the oscillations disappear, a phenomenon 
not seen in the study of the one- dimensional system. 

In the one- dimensional model the structure of the slip zone in an earthquake 
is always trivial since it is connected. In two dimensions the situation is quite 
different. The structure of the slip zone is intimately related to the structure of the 
displacement field In Figs. 10 - 12 we indicate by cross those blocks which are 
subject to an eleastic force greater than 0.9 at a particular time. The value 0.9 is 
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chosen solely for illustrative purposes. We see that the softest system is quite regular 
and the boundary between the set of blocks which experience a force in excess of 0.9 
and the rest of the blocks consistes mostly of straight lines. As we increase A the 
system becomes more irregular with valleys of all sizes. The system becomes more 
homogeneous as A is increased still further. The shape of the slip zone for typical 
earthquakes occuring in systems with these three values of the stiffness parameter 
are illustrated in Figs. 13 - 15. In the soft system the slip zone has usually a regular 
boundary, mirroring the structure of the underlying displacement field. The slip 
zone on Fig. 13 is atypical in the sense that most large earthquakes for A = 0.67 
extend over a large fraction of the system. For intermediate values of A the slip 
zone becomes more irregular. For the very stiff system, large events tend to extend 
over most of the system and the slip zone has an irregular boundary. The Hausdorff 
dimension corresponding to the slip zones in these three cases is clearly smallest for 
the most chaotic system corresponding to the interemediate value A = 0.82 for the 
stiffness parameter. 

A further qualitative insight into the dynamics of the system is obtained by 
following the development of earthquakes in real time and see it move block by 
block. In the soft system the earthquake front propagates in the beginning along 
ribbons which make a 45° angle with the coordinate axes. This is expected because 
in this way two blocks can cooperate to topple a block which they have as a common 
nearest neighbour. As the slip zone grows in size it often becomes roughly hexagonal 
for the largest events. In the stiffest system we also observe the hexagonal expansion 
of the slip zone but the ribbons are absent and earthquakes engulf the more stuck 
regions. The intermediate system is most irregular as expected. 

The average Hausdorff dimension is not the relevant measure of the structure 
of the slip zone. This is illustrated in Figs. 16 - 17 for the system with A = 0.67. 
Fig. 16 is a histogram of the Hausdorff dimension of events with more than 100 
blocks. There are two peaks in the distribution at the values 5 = 1.86 and 5 = 1.52. 
If we look at Fig. 17, which shows the distribution of events with more than 1000 
blocks participating, we see that the peak at the lower value of 5 has disappeared. 
The earthquakes in the A = 0.67 system therefore fall into two classes: Those that 
extend over almost the entire system and have a slip zone with straight edges and 
those that are bounded in size. These two classes of events correspond to the tail 
and main body of the R — fi distribution in Fig. 3. The slip zone illustrated in Fig. 
13 corresponds to the events that have a Hausdorff dimension close to 2, but most 
events in this class are in fact much larger. 

We observe a similar distribution of the Hausdorff dimension in the stiff system 
A = 0.95 with large Hausdorff dimension for the large events whereas the events in 
the intermediate system (A = 0.82) have a Hausdorff dimension that peaks around 
5 = 1.6 and does not depend sensitively on the size cutoff, see Figs. 18 - 21. 



5 The OFC model 



The model studied by by Olami, Feder and Christensen ||, see also (15|, [16], [19 
is the one we have discussed above with the force relaxation function 6 = 0. For 
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describing real earthquake faults this is probably a highly unrealistic choice of force 
relaxation since the typical stress drop in earthquakes is believed to be of the order 
of 10% or less p2| . In this model the size distribution (defined as the logarithm of 



the number of participating blocks) falls nicely on a straight line. Other notions of 
size do not yield a linear distribution. 

In order to compare the OFC model with the one we have simulated we have 
calculated the spatial as well as the temporal corelation functions, see Figs. 22 and 
23. They are qualitatively quite similar to the ones found in the Nakanishi model 
with the main difference that the autocorrelation function does not oscillate around 
zero but around a positive value, see Fig. 23. This might very well be a finite size 
effect which are known to make the analysis of the OFC model quite subtle 



6 The algorithm 

In order to simulate the temporal behaviour of an observable X we have to run the 
system for a certain time T in order to obtain resonable statistics. In most cases 
we are not interested in finite size effects and want to have the system as large as 
possible. It is therefore important how the computational time scales with N when T 
is constant. For many parameter values the average size of an event is not a function 
of TV and therefore the generation of an earthquake is an 0(1) algorithm and the 
average time between earthquakes scales like iV" 1 . To run the dynamics in an N 
block system for a time T we therefore need a number of earthquakes proportional 
to TN. If we had an 0(1) method of finding the block closest to threshold we 
could carry out the simulation in a time O(N). Using a method introduced by 



Grassberger |L9| this possible under the assumption that the distribution of the 
forces on the blocks has no delta-function peaks. Our method is quite similar. We 
use M boxes B[0], B[l], .., B[M— 1], M oc N, which can store blocks in doubly-linked 
lists. The advantage of having a doubly-linked list is that the removal of a block 
from such a list is an 0(1) operation. The indexing of the B's is understood modulo 
M. We define e = j^zi an d Qi — fi — t for i = 1, 2, .., N, where t is the physical 
time. In our implementation we also have two other boxes A\ and A 2 which can 
hold blocks. They are used to keep track of the blocks that go over threshold. We 
need two such boxes because we must keep track of the order in which the blocks 
move in an earthquake. The idea is to keep the blocks in the B boxes ordered in 
such a way that the block closest to threshold can be found in a time 0(1). 

In the simulations reported in this paper we use a global relaxation function 
and a global A, corresponding to a homogeneous system. With trivial modifications 
it is possible to use the same algorithm where each block i has its own relaxation 
function <pi satisfying |0j| < 1 and the springs between the blocks have different 
spring constants. 

In the beginning we set t = and distribute the /j's randomly — 1 < / 4 < 1. We 
put block i in -S[ L^T^J ] ^ where [^J denotes the largest integer smaller than x. The 
boxes A\ and A 2 are empty. Then we perform the following steps consecutively: 
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1. All the blocks are experiencing force between — 1 and 1, t > 0, and for % = 
1, 2, .., N block i is in B[ji] where jj = [^-\ • This implies 

-l-t<qi<l-t (21) 

and then 

A < ji < ill + M - 1 (22) 

and 

1 + t - (ji + l)e < /< < 1 + t - j t e (23) 

for all i. From this it follows that if the box -B[|_fJ] is not empty, then it 
contains the block i max closest to threshold. If it is empty, we look into the 
boxes -B[[JJ + 1]) -B[Lf J + 2],.. until we find a nonempty box, which then 
contains i max - We find the block closest to threshold and remove it from the 
box. We define f max = f imax and St — 1 — f max - We then increase t by St. This 
does not affect the ordering in the boxes for the g's will not change. After this 
f max is equal to 1 and we label the i max block i acUv e- 

2. We update the i ac tive block according to the dynamics, i.e. 

f ^active ' ^ ^P\J active) V^^/ 

and update qi acttve accordingly. Then we put it into the box B[[ - 
and remove all its neighbours i x from their boxes. Because we can adress 
them through their positions in the block-array this is an 0(1) operation. We 
modify their /'s according to 

fix ' ^ fix ^ j~y^(.f ^active ^\f iactive/) 1 (^^) 

where D is the dimension of the system and update their g's. If a neighbour 
becomes critical, i.e. has force above 1, we put it into the box A 2 . The rest is 
put back in the B boxes, i x in 

3. We look into the box Ai. If it is not empty we remove some block % from Ai, 
set iactive — i and go back to step 2. If the box A 1 is empty we go to the box 
A 2 . If A 2 is not empty we remove some block % from it and set iactive — i, 
move the rest of the blocks in A 2 to A±, and go back to step 2. If both Ai 
and A 2 are empty, then there are no more critical blocks at this time, i.e. the 
earthquake is over, and we go back to step 1. 

This completes our description of the algorithm. 

7 Discussion 

In this paper we have generalized the simulation of the Nakanishi model to two 
dimensions. We have not been able to go much beyond systems of size 200 x 200 
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so our results have considerable finite size effects which are not straightforward 
to entangle. In general these systems do not have any translationally invariant 
thermodynamic limit as we saw clearly in the case of a stiff one-dimensional system 
in jnj]. The simulations reported here also show that when the system makes the 
transition from the trivial periodic state when the blocks move independently of 
each other there is a very abrupt change with systemwide events, both in the one- 
and the two-dimensional models. 

The only qualitative difference we have found between the one- and the two- 
dimensional cases is the absence of any linear region in the R — fi distribution in 
two dimensions and the positivity of the autocorrelation function for a stiff two- 
dimensional system. 

We tried to find a transition from a one-dimensional behaviour, as one expects in 
a long and narrow system, to properly two dimensional behaviour. However, finite 
size effects did not allow us to see such a crossover clearly. In real earthquakes, the 
frequence size distribution is expected to be governed by a two-dimensional theory 
for small events but by a one-dimensional theory for large events for which the 
seismogenic zone extends through the litosphere p3| . By a careful analysis of an 



elongated two-dimensional model we should be able to see such a transition in the 
model discussed here. 
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Figure Caption 



Fig. 1. The figure illustrates, as a function of A, the average number of blocks 
participating in the 10 largest out of 2 x 10 7 events in a one dimensional system 
consisting of 10,000 blocks for values of a ranging from 2 to 5. 

Fig. 2. The figure illustrates, as a function of A, the average number of blocks 
participating in the 10 largest out of 10 8 events in a 200 x 200 block system for 
a — 3. 

Fig. 3. The figure shows the relative frequency of earthquakes as a function of 
their magnitude for A = 0.67. 

Fig. 4. The figure shows the relative frequency of earthquakes as a function of 
their magnitude for A = 0.82. 

Fig. 5. The figure shows the relative frequency of earthquakes as a function of 
their magnitude for A = 0.95. 

Fig. 6. The figure shows the relative frequency of earthquakes as a function of the 
number of participating blocks for A = 0.82. 

Fig. 7. The figure shows the relative frequency of earthquakes as a function of the 
number of participating blocks for A = 0.95. 

Fig. 8. The spatial correlation function for 3 different values of A. 

Fig. 9. The autocorrelation function for 3 different values of A. 

Fig. 10. The crosses indicate those blocks that are subject to an elastic force 
greater than 0.9 in a 200 x 200 block system with A = 0.67. 

Fig. 11. The crosses indicate those blocks that are subject to an elastic force 
greater than 0.9 in a 200 x 200 block system with A = 0.82. 

Fig. 12. The crosses indicate those blocks that are subject to an elastic force 
greater than 0.9 in a 200 x 200 block system with A = 0.95. 

Fig. 13. The crosses indicate the position of the blocks that participated in a 
typical earthquake in a 200 x 200 block system with A = 0.67. 

Fig. 14. The crosses indicate the position of the blocks that participated in a 
typical earthquake in a 200 x 200 block system with A = 0.82. 

Fig. 15. The crosses indicate the position of the blocks that participated in a 
typical earthquake in a 200 x 200 block system with A = 0.95. 

Fig. 16. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 100 blocks in a system with A = 0.67. 

Fig. 17. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 1000 blocks in a system with A = 0.67. 

Fig. 18. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 100 blocks in a system with A = 0.82. 
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Fig. 19. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 1000 blocks in a system with A = 0.82. 

Fig. 20. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 100 blocks in a system with A = 0.95. 

Fig. 21. A histogram illustrating the distribution of the Hausdorff dimension of the 
slip zone for events consisting of more than 1000 blocks in a system with A = 0.95. 

Fig. 22. The spatial correlation function in the OFC model with A = 0.8. 

Fig. 23. The autocorrelation function in the OFC model with A = 0.8. 
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